CURVATURE BLOW UP IN BIANCHI VIII AND IX VACUUM 

SPACETIMES 



HANS RINGSTROM 



Abstract. The maximal globally hyperbolic development of non-Taub-NUT 
Bianchi IX vacuum initial data and of non-NUT Bianchi VIII vacuum initial 
data is C^-inextendible. Furthermore, a curvature invariant is unbounded in 
the incomplete directions of inextendible causal geodesies. 



1. Introduction 

According to a conjecture by Belinskii, Khalatnikov and Lifshitz, the Bianchi IX 
spacetimes are a good model for the local behaviour of generic gravitational collapse, 
see Berger et al (1998) and references therein. Although they have been analyzed 
numerically, few statements concerning these spacetimes have been proven. 

There are several formulations of Einstein's equations for Bianchi IX models. One 
is the Hamiltonian description due to Misner and another is the formulation by 
Wainwright and Hsu. A brief explanation of the different variables can be found in 
Hobill et al (1994). We will use the formulation developed by Wainwright and Hsu 
(1989). In section^ we explain how these variables are obtained. 

In order to formulate the main result of this paper we need to define what is meant 
by the different types of initial data. Let G be a 3-dimensional Lie group, e^, 
i = 1,2,3 be a basis of the Lie algebra with structure constants determined by 
[ei,ej] = 7j^efe. If 7*j, = 0, then the Lie algebra and Lie group are said to be of 
class A and 

(1) 7*;- = Q.-^n'^™ 
where the symmetric matrix n'^ is given by 

(2) n^' = \iMe 



1 (i^3)kl 



Definition 1. Class A vacuum initial data for Einstein's equations consist of the 
following. A Lie group G of class A, a left invariant metric g on G and a left 
invariant symmetric covariant two-tensor k on G satisfying 

(3) Rg - hjk'^ + {tr gk)^ =0 

and 

(4) Vjtrgfc - V^hj ^ 

where V is the Levi-Civita connection of g and Rg is the corresponding scalar 
curvature, indices are raised and lowered by g. 
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Table 1. Bianchi class A. 



Type 


ni 


n2 




I 











II 


+ 








VIo 





+ 




VIIo 





+ 


+ 


VIII 




+ 


+ 


IX 


+ 


+ 


+ 



We can choose a left invariant orthonormal basis {e^} with respect to g so that the 
corresponding matrix rf^ defined in is diagonal with diagonal elements ni, n2 
and n^. By an appropriate choice of orthonormal basis ni, 712,^3 can be assumed 
to belong to one and only one of the types given in table |l|. We assign a Bianchi 
type to the initial data accordingly. This division constitutes a classification of the 
class A Lie algebras. We refer to Lemma |l| for a proof of these statements. 

Let kij = kifii, Cj). Then the matrices n*-' and fcy commute according to so that 
we may assume kij to be diagonal with diagonal elements fci , k2 and k^ , cf . (^) . 

Definition 2. Class A vacuum initial data, except type I and VIIq data with tvgk = 
0, satisfying /c2 = ^3 o-nd 77.2 = or one of the permuted conditions are said to be 
non-generic. In the Bianchi IX and VIII cases we call such data Taub-NUT and 
NUT vacuum initial data respectively. 

We will justify this definition in the remark following Lemma ^. Observe that the 
condition is independent of the choice of orthonormal basis diagonalizing n and fc, 
cf. (|^). What is meant by inextendibility is explained in the following 

Definition 3. Consider a connected Lorentz manifold (M,g). If there is a con- 
nected Lorentz manifold (M , g) of the same dimension and a map i : M —> M, 
with i{M) 7^ M, which is an isometry onto its image, then {M,g) is said to be 
C"^ -extendible and {M,g) is called a C"^ -extension of {M,g). A Lorentz manifold 
which is not C'^ -extendible is said to be -inextendible. 

Remark. There is an analogous definition of smooth extensions. Unless otherwise 
mentioned manifolds are assumed to be smooth and maps between manifolds are 
assumed to be as regular as possible. 

We are now in a position to state the main theorem. 

Theorem 1. Consider the maximal globally hyperbolic development of non- Taub- 
NUT Bianchi IX vacuum initial data or of non-NUT Bianchi VIII vacuum initial 
data. It is C"^ -inextendible and the Kretschmann scalar 

is unbounded in the incomplete directions of inextendible causal geodesies. 

Remark. One can time orient the Bianchi VIII development mentioned in The- 
orem 1^ so that all causal geodesies are future complete and past incomplete. In 
the Bianchi IX development mentioned, all causal geodesies are future and past 
incomplete. A proof, not the first, is given in Lemma 0. 
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We also obtain qualitative information concerning the asymptotic behaviour of 
solutions in the causally geodesically incomplete directions. A solution which is not 
NUT nor Taub-NUT must oscillate indefinitely in the variables of Wainwright and 
Hsu as one approaches a singularity. 

In his article concerning the global dynamics of the Mixmaster model, Rendall 
(1997) proved, among other things, the corresponding theorem for Bianchi type I 
(with tvgk 0), II, VIq and VIIq (tr^A: 0) vacuum initial data, the exceptional 
cases being the non-generic data. Bianchi type I and VIIq vacuum initial data 
satisfying tTgk — yield geodesically complete, and consequently inextendible, 
maximal globally hyperbolic developments, see Lemma |[ Thus we have 

Corollary 1. Strong Cosmic Censorship holds for Bianchi class A vacuum initial 
data. That is, the maximal globally hyperbolic development of generic vacuum initial 
data of class A is -inextendible. 

For a related paper see Chrusciel and Rendall (1995). They isolate the locally 
homogeneous vacuum initial data on compact manifolds whose maximal globally 
hyperbolic developments have smooth extensions. In the Bianchi IX case, they 
prove that the homogeneous vacuum initial data that lead to an extendible maximal 
globally hyperbolic development are the Taub-NUT initial data. They also show, 
in the compact, locally homogeneous case, that extendibility implies additional 
symmetry of the solution in the sense that the local Killing algebra has dimension 
at least four. Thus they are able to conclude that extendibility is non-generic. 

Our arguments can be extended to a more general situation. Consider initial data 
as in Definition |l| with G simply connected. If F is a subgroup of G that acts freely 
and properly discontinuously on G on the left, then we get well defined initial 
data on the quotient G/T and our results apply to them. However, one does not 
necessarily obtain all quotients of G that admit locally homogeneous data in this 
way. One can for instance take quotients of by a subgroup of the isometry group 
(of the standard flat metric) which is not a subgroup of the group of translations. 
In certain situations the process of taking the quotient yields restrictions on the 
possible initial data. These restrictions have to be analyzed for one to be able to 
answer questions concerning cosmic censorship. We refer to Chrusciel and Rendall 
(1995) for further elaboration on this point. 

Section |2| contains a sketch of a derivation of the evolution equations. The first 
form of the equations presented there is due to Ellis and MacCallum (1969). We 
then reformulate them as in the article by Wainwright and Hsu (1989) and state 
some curvature expressions we will need later. The remaining sections contain the 
proof. 

2. The equations of Wainwright and Hsu 

In this section we consider a special class of spatially homogeneous four dimensional 
vacuum spacetimes of the form 

(5) {IxG,-dt^ + x^Jit)e'»e) 

where / is an open interval, G is a Lie group of class A, Xij is a smooth positive 
definite matrix and the are the duals of a left invariant basis on G. Below, Latin 
indices will be raised and lowered by 6ij . 
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Lemma 1. Ta6Ze|^ constitutes a classification of the class A Lie algebras. Consider 
an arbitrary basis {ci} of the Lie algebra. Then by applying an orthogonal matrix 
to it, we can construct a basis {e[} such that the corresponding n' defined by 
has diagonal elements of one of the types given in table |. 

Proof. Let be a basis for the Lie algebra and n be defined as in (||). If wc change 
the basis according to e- — {A^-^)/ej then n transforms to 

(6) n' ^ {detA)-^A*nA 

Since n is symmetric we assume from here on that the basis is such that it is 
diagonal. The matrix A = diag(l 1 — 1) changes the sign of n. A suitable 
orthogonal matrix performs even permutations of the diagonal. The number of 
non-zero elements on the diagonal is invariant under transformations taking one 
diagonal matrix to another. If A = [aij ) and the diagonal matrix n' is constructed 
as in (^) we have n'^.^, = (det A)~^ Y^^i=i ^ifc*^" so that if all the diagonal elements of 
n have the same sign the same is true for n' . The statements of the lemma follow. 
□ 

Let {M,g) = {I x G, —dt^ + Xij{t)C ^ be four dimensional and G be of class 
A. Let cq — dt and Ci = a/ Zj, i=l,2,3, be an orthonormal basis, where a is a 
C°° matrix valued function of t and the Zi are the duals of Let the matrix 
A satisfy eo(^) + AB = 0, A{0) — Id where Bij —< Weoei,ej > and Id is the 
3x3 identity matrix. Then A is smooth and iS'0(3) valued and if — A^-'ej then 
< Ve(,e-,e^ 0. Assume < ^eo^ii^j >— 0. Observe that [Zi,eo] = 0. The Ci 
span the tangent space of G and < [eo, e^], cq >= 0. Let 

(7) e{X,Y)=<Vxeo,Y>, 

ep) and [ep,e^] — j^^ea where Greek indices run from to 3. The 
objects 9af3 and will be viewed as smooth functions from / to some suitable 
K*^. We get 6*00 = ^'oi = and 9ai3 symmetric. We also have 7°^ — loi — and 
7qj ~ —Oij. We let n be defined as in (||) and 

where we by abuse of notation have written tv(9) as 9. 

We compute the Einstein tensor in terms of n, a and 9. The Jacobi identities for 
Ca yield 

(8) eo(ny) - 2n^^^ai^ + ^9n,j = 0. 
The Oi-components of the Einstein equations are 

(9) cr^nfcj - n.Vfcj = 0. 

Letting bij — 2ni^nkj — tr{n)nij and sij — bij — ^ti-{b)5ij the trace free part of the 
ij equations are 

(10) eo(o-ij) + 9aij + Sij = 0. 
The fact that Rqq = yields the Raychaudhuri equation 

(11) eo(0) + %0'^=O 
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and using this together with the trace of the ij-equations yields a constraint 
(12) + (n^n'-'" - ^tr{nf) - 

Equations (||)-(|l^ are special cases of equations given in Ellis and MacCallum 
(1969). At a point to we may diagonalize n and a simultaneously since they com- 
mute (^). Rotating by the corresponding element of S0{3) yields upon going 
through the definitions that the new n and a are diagonal at to- Collect the off- 
diagonal terms of n and a in one vector v. By (|) and (|l^ there is a time dependent 
matrix C such that v = Cv so that v{t) ~ Q\/t since v{to) — 0. Since the rotation 
was time independent < Ve„ej;, Cj >= holds in the new basis. 

In order to prove curvature blow up, we need to relate the maximal existence 
interval for solutions to (||)-(|l2|) to the maximal globally hyperbolic development. 
One ingredient is the following lemma, the proof of which can be found in the 
appendix. 

Lemma 2. Let {G,g, k) he class A vacuum initial data. We can then choose a left 
invariant orthonormal basis {e[} with respect to g so that the corresponding n' given 
by is of one of the types given in table [J. We can also assume kij = k(e[, e'j) to 
be diagonal. Then there is a manifold as in where I = (t_,t+) is the maximal 
existence interval for a solution of (^-(^^ with initial data n{to) ~ n' , aij[tf)) — 
kij — tv{k)5ij /i and 9(to) = tv(k), solving Einstein's vacuum equations, such that 
the metric restricted to Mt^ — {to} x G is g and k is the second fundamental form 
of Mfg . The development has the following properties: 

1. Each My — {v} X G , V Cz I , is a Cauchy surface. 

2. // the initial data is not of type IX, the development can be time oriented so 
that it is future causally geodesically complete and past causally geodesically 
incomplete (unless the data is of type I or VIIq with tvgk = in which case 
the development is causally geodesically complete). We assume this time ori- 
entation when we speak of such developments, and always that dt is future 
oriented. The development of Bianchi IX initial data is both future and past 
causally geodesically incomplete. 

3. The causally geodesically complete developments have 9 = and atj = for 
the entire solution. The developments with one complete and one incomplete 
direction have 6 > for the entire solution. Bianchi IX developments have 
6 > in (t_,to) cind 6 < in {tQ,t^) for some to G (t_,i-|_). 

4. If the Kretschmann scalar k is unbounded as t ~f t^ then it is unbounded along 
past inextendible causal geodesies. Similarly for t ^ t^ and future oriented 
causal geodesies. 

Remark. By inspecting (^ and (|lO|) we see that if k2 = fcs and n'2 = rig then cj2 = 0-3 
and n2 = for the entire solution. Since such solutions are of NUT and Taub-NUT 
type for the Bianchi VIII and IX cases respectively, see Ellis and MacCallum (1969), 
NUT and Taub-NUT vacuum initial data yield NUT and Taub-NUT developments 
respectively. 
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Introduce, as in Wainwright and Hsu (1989), 

Stj = Bij — -B\6ij 
and define a new time coordinate r, independent of time orientation, by 

For Bianchi IX developments as in Lemma ^ we only consider the part of spacetime 
where 9 is strictly positive or strictly negative. Let = |(S22 + S33) and I]_ = 
\/3(E22 — 5]33)/2. If we let Ni be the diagonal elements of A^^^^, equations (||) and 
LOl) turn into 



N[ = (<7-4I]+)iVi 

A^^ = (g + 2E+ + 2V3I]_)A^2 

(14) = (g + 2S+ -2V3I]_)Af3 

= -(2-g)S]+-35+ 

!]'„ = -(2- g)s]_ - as*- 

where the prime denotes derivative with respect to r and 
q = 2(S^ + S^) 

(15) 5+ = i[(iV2-7V3)2_7Vi(27Vi-iV2-iV3 

5_ = ^(Ar3-_7V2)(iVi-iV2-7V3). 



The vacuum constraint (^^ becomes 

3 

4 



(16) + + -[7V2 + + - 2(iViiV2 + iV2iV3 + N3N1)] = 1. 



The equations (|l4|)-(|l6|) have certain symmetries described in Wainwright and Hsu 
(1989). By permuting Ni,N2,N3 arbitrarily we get new solutions if we at the 
same time carry out appropriate combinations of rotations by integer multiples of 
27r/3 and reflections in the (S+, S_)-plane. Below we refer to rotations by integer 
multiples of 27r/3 as rotations. Changing the sign of all the Ni at the same time 
does not change the equations. Classify points (A^i, 7V2, Af3, S+, S_) according to 
the values of A^i,-/V2,A^3 in the same way as in table ^ Since the sets Ni > 0, 
Ni < and N^ = are invariant under the flow of the equation we may classify 
solutions to (p^-(p^ accordingly. When we speak of Bianchi IX solutions we will 
assume TV^ > 0, and for Bianchi VIII solutions we will assume two Ni > and one 
< 0. 

There are three special points in the (5]+, E_)-plane. They are (—1, 0) and (1/2, ±-\/3/2) 
on the Kasner circle q = 2. They are rotated into one another by applying the sym- 
metries and the corresponding type I solutions represent isometric flat universes. If 
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the spatial topology is , the metric of these solutions can be written as 

3 

(17) ds^ = -dt^ + V t'^P^dx' (g) dx' 



i=l 



on the manifold M = M+ x R^, where M+ = e M : t > 0}, with two of the pi 
zero and one equal to 1. The rotational symmetries in the E_|_I]_-plane correspond 
to permutational symmetry of the pi. We see that if for instance p2 = Ps = we 
have a rotational symmetry in the z^x^-plane at each point of M . 

The set I]_ = 0, iV2 = ^3 is invariant under the flow of (p^)-(p^). Applying the 
symmetries to this set we get new invariant sets. These sets correspond to the non- 
generic data in Definition y. Specifically, in the Bianchi IX case they correspond to 
the Taub-NUT solutions. 



The Raychaudhuri equation (11) takes the form 

(18) 9' = -{l + q)6. 

According to Lemma || the existence interval for solutions of (p^)-( p^) is past infi- 
nite. As T — > —oo, 6{t) goes to infinity exponentially according to (|18|). 

The normalized Kretschmann scalar is given by 

(19) R = Raf5^sR"'^^^/e\ 
We have 

where, by Wainwright and Ellis (1997), the normalized E and H fields are given by 

1 2 

(20) Ei-j = - (E^ Yjkj - 2^ + 

(21) = -3S]^(,7V,)fe + TV'^'Sm*,, + ^A^'^feS,, 

and Y? — YjijY?^ /2. We may consider all matrices involved to be 3-vectors. If we 
write Eij as E = {Ei E2 E3) and similarly for H we have k = 8{\E\'^ — \H\'^). Using 
the definitions and the constraint (|l^) we compute 

AriS+ + -^(A^2-iV3)S]_ 

-^N2iY+ + a/3S_) + i(A^3 - ~ i=E_) 

^ S_(l-2S+) + (7V2-Af3)(A^2 + ^3-^i) 



(22) 


Hi 


(23) 


H2 


(24) E2 


-Es 


(25) E2 


+ E3 



3V3 

2 10 2 8 

: -S+ + + _ _ + (iV2 - A^3)2 - iVi(A^2 + A^3)- 

The fact that E and H are traceless yields the remaining components. We end this 
section with a technical lemma relating the existence intervals for the variables of 
Wainwright and Hsu and those of Ellis and MacCallum. 



Lemma 3. If for all non-Tauh-NUT Bianchi IX initial data for i\l^-^lq) we can 
prove the existence of a sequence — » —00 such that ii{Tk) does not converge to 
zero, then for every development as in Lemma^ of non-Tauh-NUT initial data, the 
Kretschmann scalar is unbounded as t t- and t —t t^ (here t+ and t- are the 
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same as the ones occurring in Lemma Similarly, if for all non-NUT Bianchi 
VIII initial data of ([7^-(|7^ we can prove the existence of a sequence — ^ — oo 
such that k{Tk) does not converge to zero, then for every development as in Lemma 
1^ of non-NUT initial data, the Kretschmann scalar is unbounded as t t^ . 

Remark. That the existence interval for solutions to (|lj)-(|l^) is past infinite is 
proven in Lemma |^. 



3. Elementary properties of solutions 



We begin by giving the past existence intervals for solutions to (|lj)-(16). Let 
g e C°°(M",M") and consider x = g{x), x{Q) — xq. Let (i-,0] be the maximal past 
existence interval for the solution. 

Lemma 4. Ift- > —oo there is no sequence tk ^ t^ , tk (^--,0], such that x{tk) 
converges. 

Proof. Assume tk t^, x{tk) X-. There is an e > and a 5 > such that we 
have a smooth flow 

$ : B,{x-) X {-5,6) 

where B^{x-) is the open ball with radius e and center x-. Choose a k such 
that x{tk) G B^{x-) and \tk - t^\ < S/2. Define y{t) = x{t) for t G {tk,0] and 
y{tk — t) = ^{x{tk),t) for t G [0, S). Then y extends x beyond the maximal existence 
interval. □ 



Lemma 5. The past existence interval for solutions to is (— oo,0]. 

Proof. For all types except IX the constraint implies S^+S'i < 1 so that Ni, N2, N3 
do not grow faster than exponentially, by (|l^). Let (t_,t+) be the maximal exis- 
tence interval for a specific solution, not of type IX. If r_ > —00 we may for any 
sequence Tk > converging to r_ extract a subsequence such that the solution 
converges, contradicting Lemma ^. In other words we may not have t_ > —00, and 
in the same way we have r+ = 00. Consider solutions of type IX. Let r < 0. We 
get 

iVi(T) = exp{ / [q{s) - 4S+(s)]ds}iVi(0) < e-^^Ni{0) 



Jo 

and similarly for N2 and N3. Thus, the Ni do not grow faster than exponentially 
as we go backward in time and by the constraint the same holds for The 
past existence interval must be (— cx3,0]. □ 

We will need a few more lemmas concerning the behaviour of solutions of type VIII 
and IX. 

Lemma 6. Consider a solution of type IX. The image 00, 0]) is con- 

tained in a compact set whose size depends on the initial data. Further, if at a point 
N3>N2> Ni and N3 > 2, then N2 > N3/IO. 



Remark. The importance of the second part of the lemma is to be found in the 
consequence that one of the Ni may not go to infinity alone. 
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Proof. Let us prove the second statement first. Assume A^3 > 2 and N3 > ION2 > 
lOiVi, so that ~N2 > -Ns/lO and similarly for A^i. We get 

(iVi - iVa)^ + ~ ^NsiNi + N2) > ^ 

contradicting the constraint (^6|). We are done. Since (TVi A'2A'3)' = 3qNiN2N3 the 
product decreases as time decreases, and we have 

(26) iNiN2N3){T) <C <oo 

for r < 0, where C is a constant. There is a compact set Ki such that if A^^ < 2, 
i = 1,2,3, then € Ki. If one of the Ni is greater than 2, we may by 

a permutation assume that N3 > N2 > Ni and by the first part of the lemma 
conclude that N2 > N3/IO. Observe that then Ni{N2 + N3) < IOC so that 

by the constraint. □ 



E^ + Ei < 1 + 15C 



Lemma 7. Consider a solution of type VIII. Assume N3 > N2 > > A^i. Then 
Nf < 4/3. If further A^g > 4 then N2 > N3/2. 

Proof. The constraint says 

(27) E^ + E^ + ^[Nf + {N2 - - 2Ni{N2 + N3)] = 1. 

Since all terms on the left are non-negative, the first part of the lemma is immediate. 
The second part follows from 

2 1 

A^a - A^a < ^ < 2 < -N3. □ 
V3 - - 2 ^ 

We end this section with a lemma needed in the sequel. Consider a positive function 
A^ satisfying A^' — hN where /i(t) — > a as t ^ —00 

Lemma 8. Assume N and h are as above. Then for all e > there is a T such 
that T < T implies 

exp[(Q; + e)T] < N{t) < exp[(Q! — €)t]. 
4. Limit characterization of special solutions 



Proposition 1. A solution to ^l^J-tjldj) satisfies 

lim (E+(T),E_(r))-(-l,0) 

r — ^ — 00 

only if it is contained in the invariant set E_ — and N2 — N3. 

Remark. The implication is in fact an equivalence. The analogue for (E+, E_) 
(l/2,±V5/2) is true due to the symmetries. 

Assume in this section that (E+, E_) converges to (—1, 0). Consider 

Observe that / is either identically zero or always strictly positive due to the fact 
that N2 = A3, E_ = is an invariant set. A function related to this one occurs in 
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Wainwright and Hsu (1989). We prove an estimate of the form /(T) < g{T,T) for 
some T and t <T and then that gir, T) — s- as r — s- — oo. Compute 

(28) /' = -(2 - q)^j:l + 4V3iVi(iV2 - Af3)S_ + 2{q + 2S+)(7V2 - . 

Lemma 9. For all e > there is a T such that 

/(T)</(T)exp[e(r-r)] 

for all T < T. 

Proof. By Lemma |^ Ni goes to zero as r — > — oo. The same is true of 2 — g and 
q + 2S+. Thus, for all e > there exists a T such that t < T implies /' < ef by 
(|8|). The lemma follows. □ 

Lemma 10. // there is a sequence — > — oo as k ^ oo such that E_|_(Tfe) < — 1, 
then S- = and N2 = N3. 

Proof. The constraint yields, in the points r^, 

J/ = S'- + |(^2 - Nsf < ^N,{N2 + m). 

Applying Lemma ^ to N1N2 and N1N3, we have, for k large enough, /(t^) < 
exp(2Tfc). Choose a finite T corresponding to e = 1 in the previous lemma. We get 
f{T) < exp(rfc + T). Letting fc ^ 00 we get f{T) — 0, but then / is identically 
zero. □ 

Lemma 11. If there is an S such that 5]_|_(t) > —1 for all t < S then S_ — and 

Proof. Using the constraint to express 2 — q in the Ni we get, if we assume — 1 < 
S+ <0, 

(S+ + 1)' < 6N^ + 6\Ni{N2 + N3)\. 

Let T be such that the right hand side < 2e'^^ and -1 < i;+(r) < for all t <T. 
Integrate the inequality to obtain 

< S+(r) + 1 < e^^ 

for all T <T. But then the constraint yields 

If - -In^ + lN,iN2 + N3) + (1 - E+)(l + S+). 

By the above argument we have control over the last term, and as in the previous 
lemma we have control over the first two terms. Thus, there exists an S' such that 
T < S' implies /(r) < e"^. We deduce, using e = ^ in Lemma |9[ that / is identically 
zero. □ 

5. Properties of past limit points 

If there is a sequence — > — 00 such that a solution to (p^-(p^ evaluated at 
converges to x, then x is said to be an a-limit point of the solution. The a-limit set 
of a solution is the union of its a-limit points. Consider a solution of type VIII or 
IX. Assume it has an a-limit point (ni, n2, 713, s+, s_). The objective of this section 
is to prove that if (s+, s_) is not one of the special points it has an a-limit point 
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at which the normahzed Kretschmann scalar is non-zero. The foUowing lemma will 
be needed, 

Lemma 12 (Kendall (1997)). Let U be an open subset ofMJ", f E C°°(]R",R") 
and let F : U ^ be a continuous function such that F{x{t)) is strictly monotone 
for any solution x{t) of 

(29) x'{t) = f{x{t)) 

as long as x{t) is in U . Then no solution of ( jj^ j whose image is contained in U 
has an a-limit point in U . 

As noted in RendaU (1997), 

Lemma 13. If a solution of type VIII or IX has an a-limit point {ni,n2-,n3^ s+, S-), 
then one of the Ui must be zero. 

Proof. We have 

liViiVaiVal' = 3(?|iViA^2A^3|. 

If g = then + > and consequently \NiN2N^\ is strictly decreasing as 
we go backward in time if it is non-zero. The statement follows by Lemma ^[ □ 

Corollary 2. If a solution of type VIII or IX has an a-limit point (ni, ?i2, 713, s+, s_), 
then s\ + s'^_ <l. 

Observe that points of type VIIq defined by E+ = -1, I]_ = 0, iV2 = ^ 
and A^i = are fixed points to ([l^)-(|l^). The same is true of the rotated points. 
These points are called flat points of type VIIq. These solutions correspond to the 
special points in the following sense. The vacuum spacetime given by the manifold 



M — R+ X W' with a metric as in (17) where pi = 1 and P2 ^ ~ ^ can be 
viewed as a type VIIq spacetime by giving a different group structure. Then 
the Wainwright and Hsu variables corresponding to a suitable frame (left invariant 
under the new operation on R'^) will be a flat point of type VIIq. 



Theorem 2 (RendaU (1997)). The a-limit set of a solution of (|_?4/)-(|_?6[) of type I, 
II, VIq or VIIq is a single point of type I or a flat point of type VIIq, the latter 
only being possible if the solution is time independent. 



Lemma 14. Consider a solution to of type VIII or IX. If it has an a- 

limit point (rii, 712, na, s+, s_) such that (s+,s_) is not special, it has at least two 
a-limit points on the Kasner circle at least one of which is non-flat. In consequence 
there is an a-limit point such that the normalized Kretschmann scalar evaluated at 
it is non-zero. 



Proof. The solution of (|14|)-(|16[) with initial value (ni, 7^2, 713, s+, s_) consists of 
a-limit points to the original solution. Let us denote it (A^i, N2, N3, S+, S-). The 
different solutions will be referred to as the solution and the original solution. Let us 
first prove that there is a non-flat a-limit point on the Kasner circle. Since (s+, s_) 
is not special we know that the solution is not a flat point of type VIIq and we can 
consequently assume that the a-limit set of the solution is a single point of type I 
by Theorem |[ say (0, 0, 0, cr+, cr_). If (cr+, (t_) is not special we are done so assume 
it is (—1,0). By arguments given in the proof of Theorem ||, Ni,N2,N3 have to 
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belong to a compact set as we go backward in time. If S_) does not converge 
to (—1, 0) we may thus construct an a-limit point to the solution different from the 
only one that exists. By Proposition]^ we conclude that E_ = and N2 = N3. If 
A^i = we get E+ = ±1 by the constraint so that = — 1 contradicting the fact 
that (s+, s_) is not special. Thus Ni ^ and N2 = = 0. Applying the flow to 
this point we conclude that limT-^oo(S+(T), S_(r)) = (1,0), again the details are 
to be found in the proof of Theorem ||. Since the a-limit set is closed we conclude 
that the original solution has a non-flat a-limit point on the Kasner circle. 

We now prove that there are at least two a-limit points on the Kasner circle. Given 
the non-flat limit point (cr+,cr_), there is a neighbourhood i?e((T+,(T-) of it in the 
I]+I]_-plane such that two Ni decrease exponentially as exp(ar) and one increases 
as exp(— ar) for some a > if (S_|_,S_) is in that neighbourhood. The solution 
must thus leave the neighbourhood as we go backward in time since otherwise 
Lemma |^ or Lemma ^ would yield a contradiction. If the solution evaluated at rj, 
converges to the non-flat limit point on the Kasner circle, we can thus construct a 
sequence Sk such that the S-|_I]_ -variables of the solution evaluated at it converges 
to a point on the boundary of {a+ , cr_ ) . We can also assume that the A^j-variables 
of the solution, evaluated at Sk, converge. The limit point will have at most one 
A^i-variable non-zero by the construction. If all the iVi-variables are zero we have 
obtained a second a-limit point on the Kasner circle. If one is non-zero, we may 
apply the flow to it to obtain two a-limit points on the Kasner circle, cf. the proof 
of Theorem ^ Since the normalized Kretschmann scalar is non-zero on a non-flat 
limit point on the Kasner circle, see (|2^), the last statement follows. □ 



6. Existence of limit points 

In order to prove that an a-limit point always exists we need a few lemmas. Let 
||7V|| denote the Euclidean norm of = (iVi, A2, iVg). 

Lemma 15. Consider a Bianchi VIII or IX solution. If \\N\\ 00 as t ^ —00 
we can, by applying the symmetries to the solution, assume that iV2, A3 00 and 
Ni,Ni{N2 + N3)^0 as T ^ -00. 

Proof. We have 

(30) IA1A2A3I <C Vr < 
where C > 1. Let tq < be such that 

(31) A^i^ _^ A^l + A| > 300C^/^ Vt < TO. 

If for r < TojAil = IA2I < IA3I we get A^3 > lOC^/^ and |A^i| = |A^2| > C^/^ by 
inequality ( |3T| ) and Lemma ^ or ^. Since this is not reconcilable with inequality ( |30|) 
we conclude that one Ni, say A^i, must satisfy jA^ij < A^2 and |Ai| < A3 Vr < tq. 
For a Bianchi VIII solution A^i < 0. But then A^| + A| < 1G1A^| by Lemma | or 
and similarly for N3. Consequently A^2, A3 — > 00 as t —00. The rest of the 
conclusions of the lemma follow by inequality (|30|). □ 

Lemma 16. Consider a Bianchi IX solution. If —1 < S]-|_(ti) < 1/3 and 9Ni < 
A2 + A^3 in [t2,ti] then S+(ti) < S+(t) for t e [r2,ri]. 
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Proof. Using the constraint ( [T6| ) we have 

(32) K = i<l- 2)(S+ + 1) + liN^ - N,{N2 + N,)). 

The constraint also yields <? - 2 < 3iVi(iV2 + iVg) so that, if -1 < E+ < 1/3, 

< ^{9N, - N2 ^ N3)N,. 

The lemma follows. □ 

Lemma 17. Consider a Bianchi VIII solution to ^l^D-^ldj) and let I = [t, s], where 
s < and belongs to the existence interval. There is an eo > with the following 
property: for all eo > e > there is an A > such that if N2{t), A^3(t) > A for all 
T £ I, then — l + 2e < S+(s) implies — 1 + e < S+(t) for all t £ I. A only depends 
on e and the solution. 

Proof. Let e > 0. We have |iViiV2iV3| < C in / where C only depends on the 
solution. Consider 

(33) ^'+ = {q- 2)(S+ + 1) + ^-{Nl - N^(N^ + N^)). 

Observe that the last term is the only one pushing in the negative direction, 
since q < 2 and A^i < for Bianchi VIII. In order for E+ to attain — 1 + e in / it 
must go between — 1 + 2e and — 1 + e in /. Due to (|l^ ) the last term of the right 
hand side of (|33| ) will then behave as exp(2T — 2Ti)r/, if e is small enough, where ti 
corresponds to — 1 + 2e and 77 may be chosen arbitrarily small by choosing A large 
enough. In consequence it is impossible for to attain — 1 + e if we choose A big 
enough depending only on e and the solution. □ 

Proposition 2. A Bianchi IX solution has an a-limit point. 

Proof. Assume ||A^|| — > 00. By Lemma |5| we may assume Ni,Ni{N2 + N3) 
and N2, Ng, 00 as r — > — cx). Since N2N3 — s- cx) we conclude that 

(34) [ (E^ + + E+)ds ^ 00 

Jo 

but then 

(35) / E+ds — !• —00. 



Since 9A^i < N2 + N3 for t < T Lemma and (|35|) yield the conclusion that 
E+ < Vt < T. 

Assume there is a T2 < T such that — 1 < E+(t2). Then S+(t2) < E_|_(t) Vt < T2 
by Lemma |l|. Let A = {t : E^ + S^^ < and At = An [r,T2]. Since 

E'^ < Vr < T2 , cf. the proof of Lemma |l^, we have by ( |3^ ) 

E+(r2) - E+(t) < / ^'^ds< (l + E+(r2)) / (q - 2)ds 

JA^ JAr 

< 2(l + E+(r2)) / (E^ + + E+)ds. 

JA-r 

By (^) the expression on the right tends to — 00 as r — > — cx) so that E+(t) 00 
which is a contradiction. 
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Assume there is an S such that E+(t) < — 1 Vt < S*. Then the constraint yields 
(—1,0) and Proposition |l| yields E_ = 0, so that all the Ni must 
decrease for r < as we go backward in time by ( p^ and (|l^). We have a 
contradiction to our assumption. 

In other words there is a sequence such that Ni{Tk) is bounded. Since q{Tk) is 
contained in a compact set by Lemma ^ we may extract a convergent subsequence 
and the proposition follows. □ 

Corollary 3. Consider a Bianchi IX solution. For all e > there is aT such that 
T < T implies 

q{T) < 2 + e. 

Remark. For a Taub-NUT solution there is a T such that (/(t) > 2 for all t <T. 

Proof. Since there is an a-limit point, Lemma |l^ proves that there is a sequence 
Tk such that {NiN2N3){Tk) but since N1N2N3 is monotone we may conclude 
that the product converges to zero. Let 10^'^ > 77 > 0. Let T be such that 
{NiN2N3){t) < rp for aU t < T and assume for a t < T that iVg > A^2 > ^i- We 
get Ni{t) < T]. If, in r, iVg > 2 Lemma | yields Ni{N2 + N3) < ?/. If TVg < 2 in r 
we get Ni{N2 + N3) < Ari. Thus we have 

g < 2 + 1277 

and the statement follows by the constraint (|l6|). □ 
Proposition 3. A Bianchi VIII solution has an a-limit point. 

Proof. Assume ||iV|| — > 00 as r— > —00. The consequences given in Lemma [l5| will 
be used freely below. As N2N^ ^ 00 we get 

(36) / + + E+)ds ^ -00 



by inspecting (|l|). If E+ -1 Proposition yields N2 = N3 and S- = 0. By the 
constraint and Lemma |[ 1 + decays exponentially so that all the Ni converge to 
finite values, contradicting our assumption. If not there is an e > and by Lemma 
|l^ a T such that t <T implies —1 + e/2 < S+(t). We get, sooner or later, 

= (9 - 2)(E+ + 1) + ^{Ni - Ni{N2 + N3)) < i^l + E^ + E+ - E+ - l)e 



+ -{N^ - Ni{N2 + N3)) < (E^ + e2_ + E+)e 



so that 



E+(r2)-E+(T) <e^ '(E^ + E^^+Eh 



)ds —00 



by (1361). But that is not possible. The existence of an a-limit point follows as in 
the previous proposition. □ 



Theorem 3. // a Bianchi VIII or IX solution to /\l^-i\lq) is not contained in the 
invariant set E_ = 0, N2 = N3 or one of the sets found by applying the symmetries 
to it, the solution has an a-limit point (m, 712, ns, s+, s_) for which (s+,s_) is not 
special. 
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Proof. Consider a solution not contained in the invariant sets mentioned and let 
(ni, n2, ^3, s_|_, s_) be an a-limit point. If (s+,s_) is not special we are done, so 
assume (s+,s_) = (—1,0). Let Tk be a sequence such that the solution evaluated 
at Tk converges to the given a-limit point. Since (E+, E_) does not converge to the 
special point there is an < e < 10^^ such that for k large enough we may for each 
Tk find the first time, call it Sk < Tk, that (I]+(sfc), S_(sfc)) £ dBf{—l,0). If there 
is a subsequence of Sk such that Ni{sk) is bounded we are done so assume not. 

The goal is to go backward in time starting at Sk in order to produce a sequence 
with the desired properties. Since q is bounded as we go backward in time we must 
have Tk — Sk ^ oo so that Ni{sk) since N[ = {q — 4E+)A^i and q — 4E+ is 
roughly 6 in B,{-1, 0). We get A^3(sfc), A^2(sfc) oo and iVi(iV2 + N^) ^ in Sk. 

Bianchi IX. There is an 77 > such that —1 + rj < E_|_(sfe) < by the constraint 
if k great enough, and then we also have 9Ni{sk) < ^2(sfc) + ^3(sfc)- If the latter 
inequality holds Vr < Sk we get a contradiction to the existence of our limit point 
by Lemma ^ Let Vk < Sk be the first time 9A^i = A^2 + ^3 or one of N2, 
becomes equal to A^i. Since N1N2N3 0, cf. the proof of Corollary ||, we may 
assume 

(37) N1N2N3 < 1. 

If Niivk) = N2{vk) < Nsivk) then Ni{vk) = N2{vk) < I by and N^ivk) < 10 
by Lemma H and similarly if N2 and N3 change roles. If Ni{vk) < Ni{vk) i — 2,3 
and 9Ni{vk) = N2ivk) + N^^Vk), then Ni{vk) < 1 and N,{vk) < 9, t ^ 2,3. 
Regardless, we have Ni{vk) < 10. For k great enough there is a tk and a Uk with 
Vk < Uk <tk < Sk such that Nsitk) = 10^°, Nsiuk) = 10^" and 

(38) IQi" < iV3(r) < 10'° Vr e [uk,tk]. 
In consequence 

(39) 10^ < A^2(r) < lO'i Vr e [uk,tk] 

by Lemma 1^. The inequalities above are used to prove that tk — Uk > 1 from which 
we deduce the existence of an g [uk,tk] such that |E_(rfc)| < 1/10. That will 
suffice to prove the lemma. 

We get 

20<ln^^4H= / '(g + 2E+-2V3E_)ds. 



But Ni{N2 + A^s) < 2 • 10"" in [uk,tk] by (|37|)-(|3g|) so that g < 3 in that interval 
by the constraint (|l|), and thus tk - Uk > 1. We have \N2lN'i - 1| < 2 • 10"^° in 
[uk,tk\ by the constraint and the above. If |S_| > 1/10 in [uk,tk] we get 



But 



I / 4V3E_dr| > -. 

Jut ^ 



Mtk) ['\ /^^ ..^N2{uk) 



exp( / 4\/3E_dr)- 

•lui. 



Nz{tk) 'N^{uk) 

and we have a contradiction. Consequently there is an G [uk,tk] such that 
-1 + < S+(rfc), |E_(rfe)| < 1/10 and N,[rk) < lO^i. The conclusion of the 
theorem follows. 
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Bianchi VIII. There is an 77 > such that -1 + 2?/ < i;+(sfe) < 0. Let A = A{ri) > 1 
be as in Lemma ^ and 10 A < N3{t) Vr e [wfe, s^] and N3{vk) = 10 A, observe that 
by the existence of the hmit point the Ni must become small. Then we may argue 
as in the Bianchi IX case. □ 



7. Conclusions 

Theorem 4. Consider a Bianchi VIII or IX solution of (|I^-(|7^ which is not of 
NUT or Taub-NUT type respectively. Then there is a sequence Tk — > —00 such that 

kiTk) ^ C ^ 0. 

Proof The Taub-NUT and the NUT solutions correspond to I]_ = and N2 = N3 
or one of the sets obtained by applying the symmetries. By Theorem |^ and Lemma 
^the theorem follows. □ 

Theorem 5. For a Bianchi VIII or IX solution to J7^J-f|7^j which is not a NUT 
or Taub-NUT solution (S+,I]_) cannot converge to a special point on the Kasner 
circle and the a-limit set contains at least two distinct points of type I, at least one 
of which is non-flat. 

Proof. The first statement is Proposition |^ and the last statement follows from 
Theorem || and Lemma □ 

Observe that this theorem says that the solution does not converge; the shear 
variables will oscillate indefinitely. 

Proof of Theorem M. Let (M, g) be the globally hyperbolic Lorentz manifold ob- 
tained in Lemma 0. Assume there is a connected Lorentz manifold (M, g) of the 
same dimension and a map i : M ^ M which is an isometry onto its image with 
i{M) 7^ M . Then there is a p G A/ — i{M) and a timelike geodesic 7 : [a, h] M 
such that 7([a, 6)) C i{M) and 7(6) — p. Since 7|[a,b) can be considered to be a 
future or past inextendible timelike geodesic in M either it has infinite length or 
the Kretschmann scalar blows up along it, combining Lemma |^ and Theorem 
^. Both possibilities lead to a contradiction. Observe that this proves that (M, g) 
is the maximal globally hyperbolic development. That the Kretschmann scalar is 
unbounded in the incomplete directions of inextendible causal geodesies also follows 
from Lemma ^ || and Theorem ^ □ 



Appendix 

We here prove the technical results of Lemma || and |3[ 

Proof of Lemma |^. We begin by proving that we obtain a solution to Einstein's 
vacuum equation with the correct initial conditions. Let e^, i = 1,2,3 be a left 
invariant orthonormal basis. We can assume the corresponding n' to be of one of 
the forms given in table |l] by Lemma |l|. The content of is that kij = fc(e';, e'j) and 
n' are to commute. We may thus also assume kij to be diagonal without changing 
the earlier conditions of the construction. If we let n(to) = n', 0(io) = tr^fc and 
cTijfto) = kij — 05ij/'^ then is the same as (|l^). Let n, a and 9 satisfy (||), 
( p^ and (^1]) with initial values as specified above. Since (|l2|) is satisfied at to it 



is satisfied for all times. For reasons given in connection with (12) n and a will 
remain diagonal so that (||) will always hold. 
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Let M — I X G where / is the maximal existence interval for solutions to (g)- 
(IH). We construct a basis e^, define a metric by demanding that the basis be 
orthonormal and show that the corresponding n, a and 9 coincide with n, a and 9. 
We will thereby have constructed a Lorcntz manifold satisfying Einstein's vacuum 
equations with the correct initial conditions. 

Let Tii and cr^ denote the diagonal elements of n and a respectively. Let /i(<o) = 1 
and = 2a, - 9/3. Let = {f2hY'\ as = {hh?'\ as = (/i/s)^/' and 

define = a^e^. Then h associated to equals n. We complete the basis by 
letting eo — dt- Define a metric < • , • > on Af by demanding Cq, to be orthonormal 
with < eo,eo >= —1 and < 6^,6^ >= 1, i — 1,2,3 and let V be the associated 
Levi-Civita connection. Compute < WeoSi,ej >— 0. If 9{X,Y) —< \/xeo,Y > 
and — 9{e^, e^,), then (?oo = ^io = ^oi = 0. Furthermore, 

— eo(aj)(5ij = — 

(no summation over j) so that 6*^- is diagonal and tr6' = 9. Finally, 

The constructed Lorentz manifold thus satisfies Einstein's vacuum equations. Next 
we prove that each = {w}xGisa Cauchy surface. The metric is given by 



where are the duals of e^. A causal curve cannot intersect twice since the t- 
component of such a curve must be strictly monotone. Assume that 7 : (s_, s+) — *■ 
M is an inextendible causal curve that never intersects My. Let t : M I he 
defined by t[(s, ft,)] = s. Let sq G (^-,5+) and assume that i(7(so)) = ti < v 
and that < 7',9t >< where it is defined. Thus i(7(s)) increases with s and 
f(7([so, s+))) C [ti, v\. Since we have uniform bounds on from below and above 
on [ti , v\ and the curve is causal we get 



on that interval, with C > 0. Since 

, , r+ o 7 , 

— < 7 , eo > as = / — ; — as < v — t\ 

the curve 7|[so,s^), projected to G, will have finite length in the metric p on G defined 
by making an orthonormal basis. Since p is a left invariant metric on a Lie group 
it is complete and sets closed and bounded in the corresponding topological metric 
must be compact. Adding the above observations, we conclude that 7([so,s+)) is 
contained in a compact set. For each sequence Sk — > s+— there is thus a subsequence 
s„fc such that 7(snfc) converges. Since t(7(s)) is monotone it converges. We cannot 
have two limit points since that would contradict the causality of 7. Thus 7 must 
converge as s ^ s+ so that it is extendible. By this and similar arguments covering 
the other cases, we conclude that M„ is a Cauchy surface for each v £ (f_,t+). 
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Next we prove the statements made in Lemma ^ concerning causal geodesic com- 
pleteness. Let us first time orient the different manifolds. Consider manifolds which 
are not of type IX. Then the constraint ( p^ ) yields the inequality CTy CT*-' < 26'^/3, 
cf. (p7|). Then the Raychaudhuri equation yields |eo(^)| < (P' ■ Consequently, if Q 
is once zero it is always zero. By considering equations (||)-(12) one concludes that 
= is only possible if aij — n.ij = or if two of the diagonal elements of are 
equal and constant, the third is zero and = 0. Time orient the manifolds which 
are not of type IX and which do not have 9 ~ Q hy demanding that 6 be positive. 

Observe that for a manifold which is not of type IX, 6 decreases in magnitude with 
time, so that it is bounded to the future. By the constraint the same is true 
of aij . Using (^) we get control of and conclude that the solution may not blow 
up in finite time. The interval / = (t_, t+) in Lemma ^ must thus have = oo. 

We now prove future causal geodesic completeness for manifolds which are not of 
type IX. Let 7 : (s_,s+) ^ M be a future directed inextendible causal geodesic. 
We prove that s+ must be infinite. Let the function t be defined as in the previous 
lemma. Since every A/^,, w e / is a Cauchy surface, t(^{s)) must cover the interval 
/ as s runs through (s_,s-|_). Furthermore, i(7(s)) is monotone increasing. Define 

Let So G {s-,s+) and compute 

(40) / ~fo{u)du^t{-f[s))-t{^{so)) 

J So 

so that the right hand side goes to 00 as s — > s+. If we can prove that /o 
may not become unbounded in a finite s-interval, we are done. If = 0, then 
(t_,i+) = (—00,00) and /o is constant so that all causal geodesies are future and 
past complete. We exclude this case from now on. Compute 

k=l 

where Ok are the diagonal elements of Oij . Consider functions of t as functions of s 
by evaluating them at ^(7(5)). Compute, using Raychaudhuri's equation (ITTI), 

k=l k=l 

where are the diagonal elements of fiy . 

3 /o\V2 



E 



1 




using the tracelessness of aij. By making a division into the three cases X]fe=i — 
6*2/3, 612/3 < X;Li '^k < 26'V3 and 26*2/3 < X^Li ^fc and using the causality of 7 
we deduce 



(41) |(/„.) > '-^S'fl 

Observe that this estimate holds for a Bianchi IX solution as well. 
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Since fo9 is negative on [sq, s+), its absolute value is thus bounded on that interval. 
If S-|_ were finite, 9 would be bounded from below by a positive constant on [sq, s+), 
since 

\f\<~fo0'<C0 
as 

on that interval for some C > 0. Then (Ell) would imply the boundedness of /o on 



[so,s+). That would contradict (40) 



Excluding Bianchi IX for the moment, the inequality ( ^ ) also proves that all causal 
geodesies are past incomplete since it implies that fo0 must blow up after a finite 
s-time going into the past. 

Consider a Bianchi IX manifold. By a theorem by Lin and Wald (1989a), the trace 
of the second fundamental form will be zero for some to G I — {t-,t+) (in their 
paper they demand that G have topology S^, but this is not necessary for their 
argument). However, on J_ — (t-,to) and on /+ — (to,t-i_)6 ^ since it is only 
zero once (this can be seen by inspecting equations (||)-(^2|)). A future directed 
inextendible causal geodesic 7 in a Bianchi IX manifold must consequently have 
points Si,S2 with (^(^(7(5))) positive (negative) for s < s\ (s > 52)- The inequality 
( ^l| ) now implies that the geodesic cannot be defined more than a finite interval 
before S\ nor more than a finite interval after S2- We have future and past causal 
geodesic incompleteness. We also conclude that i_ > —00 for all manifolds that do 
not satisfy = and t+ < 00 for Bianchi IX manifolds since the curve defined by 
7(5) — (s, e) is a geodesic. 

Finally, we prove the last statement of Lemma ^. Assume the Kretschmann scalar 
is unbounded as t ^ i_ and let 7 be a past inextendible causal geodesic. Since 
each is a Cauchy surface, 7 must pass through each of them and thus the 
Kretschmann scalar must be unbounded along it. □ 

Proof of Lemma Consider a Bianchi IX solution of (p|)-(|l^. As observed in 
the proof of Lemma |^ the existence interval I = {t-, t+) can be divided into /_ = 
{t-,to), /+ = (<o, t+) and to, where to is the only zero of 9 in I. We now relate the 
different time coordinates on J_. 

According to equation (13) r has to satisfy dt/dr = 3/9. Define rit) = J^* 9{s)/3ds, 
where ti G /_. Then t : I- t{I-) is a diffeomorphism and strictly monotone on 
I- . By equation ( pT| ) 9 decreases so that it will be positive in /_ and r will increase 
with t. 

Since 9 is continuous beyond to it is clear that r(t) ^ tq £ K as i ^ to- To 
prove that t ^ t- corresponds to r — > —00 we make the following observation. 
One of the expressions 9 and d9/dt is unbounded on (t_,ii], since if both were 
bounded the same would be true of aij and rnj by ( pT| ) and (^) respectively. Then 
we would be able to extend the solution beyond i_, contradicting the fact that / is 
the maximal existence interval (observe that i- > —00 by the proof of Lemma 
If T were bounded from below on then 9 and 9' would be bounded on T((t_, ti]) 
by Lemma |5| and thus 9 and d9/dt would be bounded on (t_,ti]. Thus t t^ 
corresponds to r — > ~oo. Since — > cx) as r — > —00 by (|l^), K{Tk) — K{Tk)9^ (Tk) is 
unbounded. Consequently 

limsup\Rc,p^sR°''^^^\ = 00 
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for all non-Taub-NUT initial data given the assumptions of the Lemma. Similar 
arguments yield the same conclusion for t ^ and for < t_ given non-NUT 
Bianchi VIII initial data. □ 



Acknowledgments 



The author would like to express his gratitude to his advisor, Lars Andersson, for 
suggesting the problem and helpful discussions, and to Alan Rendall for reading 
the argument and suggesting improvements. 



References 



Berger B K. Garfinkle D . Isenberg J, Moncrief V and Weaver M 1998 Mod. Phys. Lett. A13 
1565-74, ^r-qc/9805063| 

Chrusciel P and Rendall A 1995 Ann. Phys. 242 349-85, lgr-qc/9410040 
Ellis G and MacCallum M 1969 Commun. Math. Phys. 12 108-41 

Hobill D, Burd A and Coley A 1994 Deterministic Chaos in General Relativity (New York: 
Plenum) 

Lin X-F and Wald R 1989a Phys. Rev. D 40 3280-86 

Lin X-F and Wald R 1989b Phys. Rev. D 41 2444-48 

Rendall A 1997 Class. Quantum Grav. 14 2341-56, gr-qc/970303e 
Wainwright J and Hsu L 1989 Class. Quantum Grav. 6 1409—31 

Wainwright J and Ellis G F R editors 1997 Dynamical Systems in Cosmology (Cambridge; 
Cambridge University Press) 40 



Department of Mathematics, Royal Institute of Technology, S-100 44 Stockholm, Swe- 
den 



